Ordering in the pyrochlore antiferromagnet due to Dzyaloshinsky-Moriya interactions. 
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The Heisenberg nearest neighbour antiferromagnet on the pyrochlore (3D) lattice is highly frus- 
trated and does not order at low temperature where spin-spin correlations remain short ranged. 
Dzyaloshinsky-Moriya interactions (DMI) may be present in pyrochlore compounds as is shown, 
and the consequences of such interactions on the magnetic properties are investigated through mean 
field approximation and monte carlo simulations. It is found that DMI (if present) tremendously 
change the low temperature behaviour of the system. At a temperature of the order of the DMI 
a phase transition to a long range ordered state takes place. The ordered magnetic structures are 
explicited for the different possible DMI which are introduced on the basis of symmetry arguments. 
The relevance of such a scenario for pyrochlore compounds in which an ordered magnetic structure 
is observed experimentally is dicussed. 

PACS numbers: 75.10.Hk, 75.30.Et, 75.30.Gw, 75.50.Ee 



I. INTRODUCTION 

Frustration in magnetic systems can lead to unconven- 
tional magnetic ground states and peculiar low tempera- 
ture behaviours For example, in two-dimensional and 
three-dimensional highly frustrated systems, the conven- 
tional Neel ground state may be destabilized if frustration 
is associated to a weak connectivity of the lattice. This 
is the case for example in the kagome and the pyrochlore 
lattices (see Fig. In these lattices, in the presence 
of nearest neighbour antiferromagnetic interactions, the 
ground state is a disordered spin liquid state character- 
ized by the abscence of magnetic long range order. 

The classical Heisenberg model on these two lattices 
has a macroscopic degeneracy of the ground state Q 
which prevents any Neel-like ordering at T = 0. At small 
but finite temperature, there is a tendancy to coplanar 
order in the kagome lattice Q but no order of any kind 
is observed in the pyrochlore lattice |^ . For the py- 
rochlore lattice, in the extreme quantum limit of spins 
S' = i, there is clear evidence for the absence of long 
range spin-spin correlations (and thus the abscence of 
Neel-like ordering) 0. However, the nature of the disor- 
dered state (resonating valence bond spin liquid, valence 
bond crystal,. . . ) is still unclear Isl (9|, [ij, [ill []J| . 

On the experimental side, the pyrochlore antiferromag- 
nets either do not order at low temperature, or have a 
critical temperature much smaller than the Curie- Weiss 
temperature Ocw^ which is a general trend for highly 
frustrated compounds |3 . 

The reason why a number of pyrochlore compounds 
order at low temperature is that the spin liquid state is 
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FIG. 1: The pyrochlore lattice. The interplay of the frustra- 
tion of the tetrahedral unit cell with the weak connectivity 
of the (corner sharing) tetrahedra provides peculiar magnetic 
properties to the pyrochlore antiferromagnet. 



very sensitive to any additional term in the hamiltonian 
such as on-site anisotropics , dipolar interactions [l^ , 
biquadratic interactions, next nearest neighbour interac- 
tions. . . Indeed, nearest neighbour antiferromagnetic in- 
teractions leave a high degeneracy of the ground state, or 
at least a high density of states at low energy. If the ad- 
ditional interactions are able to lift this degeneracy, they 
will be primarily responsible for the low temperature be- 
haviour of the system. 

In this article, we consider the effect of small 
Dzyaloshinsky-Moriya interactions (DMI) in addition to 
nearest neighbour antiferromagnetic interactions. The 
same study was done on the kagome lattice ^-^^ i^^ it 
will be clear from the next sections) similar consequences 
of DMI were found, although the magnetic structure is 
different in the two cases. In this paper, spins are treated 
as classical variables whereas a recent work Il6l considers 
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the extreme quantum limit of spins S=^, using another 
method. Unhke the present approach, quantum fluctu- 
ations are taken into account, but other approximations 
are made '16] . Both approaches give to some extent simi- 
lar (although not identical) results, which makes the over- 
all picture of DMI induced ordering more reliable since 
the two methods start from different limits and use dif- 
ferent approximation schemes. 

In section the microscopic origin of DMI is pre- 
sented and it is shown how to take them into account for 
pyrochlore systems, in accordance with the symmetry of 
the lattice. Section ITTll explicits the two magnetic struc- 
tures obtained for the two possible DMI we have consid- 
ered, and the possible link with experimental results is 
done. 



II. DMI AND MORIYA'S RULES FOR THE 
PYROCHLORE LATTICE 

Taking into account Coulomb repulsion and Pauli's 
principle, Anderson ^3 has explicited a microscopic 
mecanism which leads to isotropic super-exchange inter- 
actions (JS1.S2) in insulators. The main idea is that the 
energy is lowered if the spins of the electrons on neigh- 
bouring sites are in oposit directions because they can 
then minimize their kinetic energy by delocalizing on the 
nearest atoms, which is forbidden by Pauli principle if the 
spins are parallel. This fact is translated in an effective 
magnetic hamiltonian : JS1.S2. 

Taking into account a weak spin-orbit coupling (AL.S) 
and expanding in powers of A, Moriya 0, showed 
that the effective magnetic hamiltonian between two 
spins Si and S2 has the following expression : 

F„,ag = JS1.S2 + D12. (Si X S2) + Si.r.S2 + . . . (1) 

with J cx A" whereas D oc X and F which is a symmet- 
rical tensor is oc A^. These are the leading orders in A 
and the second term of is the Dzyaloshinsky- Moriya 
interaction (DMI). Despite their higher order in A, the 
third term and higher order terms should not always be 
neglected compared to the DMI for two reasons. First, 
taking all the orders into account can restore the rota- 
tionnal symmetry of the Hamiltonian as was shown in 
|20l l2lj| . This requires however some particular symme- 
try properties of the lattices which do not hold for 
the pyrochlore lattice. The second case where higher or- 
der terms should not be neglected is when the magnetic 
structure induced by the isotropic interaction (J) is essen- 
tially coUinear (ferromagnetic or antiferromagnetic) since 
the cross products Si x Sj will then be of the order of 
~ A, and the DMI term in (Q) will be of order ^ A^, the 

same order as the next term (F in eq. QJ). However, 
collinear structures do not play any particular role in the 
pyrochlore lattice |^ and we shall only consider the first 
two terms of in this paper. 



Moriya's microscopic derivation of the DMI is only 
valid for insulators but other possible microscopic mecan- 
ism relevant for other materials were explicited. for in- 
stance in systems with RKKY interactions |23, 123 . 

Whatever the microscopic origin of the DMI is, there 
are always symmetry constraints on the possible D vec- 
tors which may appear in the hamiltonian. Indeed, the 
hamiltonian must be invariant under the symmetry op- 
erations of the crystal and this will restrict the possible 
D vectors to those for which the expression is in- 
variant (under these symmetry operations). This way of 
constraining the D vectors has been given the name of 
Moriya's rules We will make use of these rules to 

determine the possible DMI in the pyrochlore lattice. 

One of Moriya's rules states that if the point between 
two magnetic sites is a center of inversion, then there 
can not exist a DMI. For that reason DMI are absent 
in crystals with a high symmetry, but this does not rule 
out DMI in the pyrochlore lattice which has no inversion 
center at the middle point between two sites. Carrying 
on the symmetry analysis one can in fact determine com- 
pletely the direction of the D vectors as follows. Consid- 
ering a single tetrahedron (the pyrochlore lattice is an 
assembly of corner sharing tetrahedra), the plane which 
contains two sites and the middle point of the opposite 
bond in the tetrahedron is a mirror plane (these planes 
are {110} planes in the cubic cell of figure 0). Applying 
Moriya's rules then implies that the D vector can only 
be perpendicular to this plane, or equivalently parallel to 
the opposite bond as shown on figure|21 There are so two 
possible DMI between two sites which correspond to the 
two directions for the D vector (and keeping the same 
order for the cross-product Si x Sj). 

The symmetry analysis (Moriya's rules) can not give 
more information about DMI. Indeed, the sign of D 
along the direction previously determined depends on the 
microscopic details of the particular compound |l3 , such 
as the number of electrons on the magnetic sites, the oc- 
cupied orbitals, the crystal field... At present, there is no 
general scheme to determine D quantatively, except in 
simple models which are unrealistic for the compounds 
we consider later on. For this reason, we will consider 
the only two possible DMI as determined by symmetry. 

One can equivalently describe the two possible DMI 
with a fixed D vector but changing the order of the cross 
product since : D.(Si x Sj) ~ — D.(Sj x S^). This is 
done without any loss of generality since changing simul- 
taneously the direction of D and the order of the cross 
product leaves the whole interaction unchanged. In the 
following the signs of the D vectors are changed while 
keeping always the same order for the cross products 
Si X Sj. These two possibilities are refered to as "di- 
rect" and "indirect" DMI and are explicited in figure [5] 

Once one particular DMI between two spins is fully 
specified, all the other DMI in the lattice are also fixed 
by symmetry. For instance the DMI within one tetrahe- 
dron can be deduced one from another thanks to the four 
three- fold rotation axes. The result is shown in figure |21 
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FIG. 2: D vectors for the DMI in the pyrochlore lattice. The 
convention is taken to fix the order for the cross products (al- 
ways D.Si X Sj with j > i). The two possible DMI are those 
obtained by varying the direction of the D vectors (D — > — D). 
The case with the D as represented on the left is refered to 
as the "direct" case and the other case (on the right) is the 
"indirect" case. Once the DMI between two spins is fully 
specified, the others DMI are also fixed and obtained by ap- 
plying the different ^ rotations around the cube's diagonals 
which leave the tetrahedron invariant. The DMI in the rest 
of the lattice are also fixed and obtained by applying appro- 
priate symmetry operations of the pyrochlore lattice (see text 
and figure I^J. 
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FIG. 3: Symetry rules to deduce the D vectors from one 
tetrahedron to its neighbours. The point where the spin Sc is 
located is a center of inversion which leaves the whole lattice 
unchanged, and the hamiltonian must be invariant in this 
symetry operation. S2 is the image of Si and thus given the 
interaction Dci.(Sc x Si), the interaction Dc2.(Sc x S2) is 
easily deduced. The two D vectors are represented on the 
figure. 

The DMI in the rest of the lattice are obtained by apply- 
ing appropriate symmetries of the lattice as follows. Two 
neighbouring tetrahedra in the pyrochlore lattice are cor- 
ner sharing and are transformed one into another by the 
inversion center located at the common site (site Sc on 
figure|3I) ■ The DMI hamiltonian must be invariant in this 
symmetry operation, and thus the DMI on one tetrahe- 
dron are easily obtained from the DMI on a neighbouring 
tetrahedron (see figure OJ. 



III. MAGNETIC PROPERTIES 

In this section, the consequences of DMI on the mag- 
netic properties of a pyrochlore antiferromagnet are stud- 



ied within mean field approximation and by classical 
monte carlo simulations. Monte carlo simulations indi- 
cate a phase transition to a long range ordered magnetic 
structure even for very weak DMI compared to antifer- 
romagnetic isotropic exchange (J). This is seen on the 
specific heat which seems to show a singularity when in- 
creasing system size (see figure^)) as well as on snapshots 
of the magnetic structure at very low temperature. Thus, 
the low temperature structure is in deep contrast with the 
spin disordered state in the abscence of DMI 0, Q . We 
start in section lTlI Al with the results in the case of direct 
DMI. The case of indirect DMI is considered in section 
ImBl The ordering of some pyrochlore compounds which 
could be due to DMI is then discussed. A slightly differ- 
ent geometry of the D together with some additional sec- 
ond nearest neighbour antiferromagnetic interactions is 
considered in relation with experimental results on min- 
eral paramelaconite (CU4O3) |[24|. 

A. direct DMI 

We start with the case of direct DMI (see sectionllllfor 
a definition of "direct" ) . We first show that DMI lead to 
a phase transition at low temperature, and then explicit 
the ordered magnetic structure. 



1. Phase transition 

Monte carlo simulations indicate a phase transition 
from a high temperature paramagnetic phase to a long 
range ordered state. This is seen in the specific heat (fig- 
ureEJ and also from snap shots of the magnetic structure 
at very low temperature. The critical temperature was 
found to be roughly of the order of magnitude of DMI, 
although it depends quite strongly on the type of DMI 
("direct" or "indirect"). 

Although the DMI are not the more intense interac- 
tions in the system, this is not a surprising result for the 
following reasons. In the presence of only nearest neigh- 
bour antiferromagnetic interactions (which are the dom- 
inant interactions), the pyrochlore lattice has a macro- 
scopic degeneracy of its ground state on a mean field level 
3] . This high degeneracy is not lifted by thermal 0, |^ 
fluctuations and the resulting low temperature magnetic 
structure is a disordered spin-liquid state. In the presence 
of DMI, even at a mean field level this macroscopic de- 
generacy is completely lifted. The DMI are thus the only 
interaction responsible for the ordering and this leads to 
the unusual behaviour that Tc ^ D even though DMI 
are an order of magnitude smaller than J (the isotropic 
exchange). This conclusion is also valid in the case of 
"indirect" DMI, although in this case there is a more 
subtle entropic effect at low temperature as explained in 
section IIII Bl Furthermore, in both cases the magnetic 
structures minimize the isotropic exchange interactions, 
so there is no competition between DMI and isotropic ex- 
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FIG. 4: Specific heat for tlie "direct" (up panei) and tlie "indi- 
rect" (down panei) cases, obtained in monte carlo simulations 
for different sizes of clusters, -j — 0.1 and N is the total num- 
ber of spins. The critical temperature and more generally the 
phase transition depends on the type of DMI (note the dif- 
ferent scales both in temperature and specific heat ranges be- 
tween the two cases). For each point with N=128 and N=1024 
(resp. N=3456), 10^ (resp. 10^) monte carlo steps were dis- 
regarded for thermalization and then 10^ (resp. 2.10®) steps 
were performed to measure the specific heat. The lines are a 
guide for the eyes. 



change. It is thus expected that Tc is mainly determined 
by DMI. 



2. Magnetic structure 

In the case of direct DMI (see sectionllllfor a definition 
of "direct"), the low temperature magnetic structure is 
unique except for the time-reversal degeneracy. All the 
spins point towards the center of the tetrahedron or all 
of them point in the opposite direction and these two 
structures are usually named "all-in" and "all-out". In 
the pyrochlore structure, the four directions of the spins 
correspond to the four [111] directions in the underlying 
cubic (fee) lattice (see figure EJ- The time reversal de- 
generacy is still present in the pyrochlore structure, but 




FIG. 5: Magnetic structure of the ground state of the py- 
rochlore lattice in the presence of "direct" DMI. The two 
structures "all-in" and "all-out" alternate from one tetrahe- 
dron to its neighbours. 



this is only a global degeneracy : the whole structure is 
uniquely determined as soon as the structure of one tetra- 
hedron is fixed to one of the "all-in" or "all-out" states. 
The resulting structure is represented on figure |S1 where 
the two structures "all-in" and "all-out" alternate from 
one tetrahedron to its neighbours. 



3. Pyrochlore compounds 

The magnetic structure obtained in the presence of di- 
rect DMI (figure was observed experimentally in the 
pyrochlore compound FeFa j25l |26| . In this compound, 
the Fe'^+ ions have a half filled electronic d shell and an 
isotropic charge distribution with zero orbital momentum 
(L=0). This is unfavorable to the appearance of sin- 
gle ion anisotropics, and in any case can not reasonably 
explain the observed phase transition at a temperature 
of 15 K. However, DMI are not the only possible ori- 
gin for this structure, third neighbour interactions 
also lead to this ground state. It is hard to discrim- 
inate between the two scenarii given the experimental 
and theoretical available predictions. However, if DMI 
arc indeed present, this should show up as an effective 
anisotropy which is not the case with Heisenberg third 
neighbour interactions. It is interesting to note that the 
parent compound NH4Fe^^Fe^^^F6 also orders at low tem- 
perature and has a magnetic structure similar to the one 
obtained with indirect DMI. This is detailed in section 

ChbI 
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B. Indirect DMI 

The consequences of indirect DMI are a bit more subtle 
to study than the direct case due to the occurence of 
entropic effects (or "order by disorder" ^^). These are 
not essential to understand the ordering of the system 
which can be predicted qualitatively on a mean field level, 
but they are important to predict the exact magnetic 
structure obtained at low temperature. 

1. Magnetic structure from monte carlo simulations 

Monte carlo simulations in the presence of indirect 
DMI indicate a phase transition to a long range mag- 
netic order at low temperature (see the specific heat on 
figure The low temperature magnetic structure is 
a coplanar structure described by a q = wave vector. 
The magnetic moments lie in either the (xy) , (zx) or (yz) 
plane, which represents a global degeneracy. The mag- 
netic structure of one tetrahedron for the (xy)-coplanar 
state is represented on figure IHI Similar structures are 
easily deduced for the (zx) and (yz)-coplanar structures. 

2. Mean-field 

At the very outset, there is no reason to solve the prob- 
lem within mean-field approximation at T = since one 
expects monte carlo simulations to be much more reli- 
able, especially in the treatment of thermal fluctuations. 
However, the merit of the following mean-field approach 
is not to reproduce (among others) the magnetic struc- 
tures found in monte carlo simulations, but precisely to 
predict magnetic structures which are not obtained in 
the simulations. This underlines that the magnetic struc- 
tures observed in the simulations are the result of an en- 
tropic selection by thermal fluctuations (which are absent 
in the mean field treatment). 

The structures found in the monte carlo simulations 
are q = structures : the magnetic elementary cell is 
identical to the crystallographic one (one tetrahedron). 
This justifies the following mean field approximation 
which assumes a wave vector q = and will give the 
corresponding magnetic structures which minimize the 
energy. The following mean field approach thus neglects 
any thermal fiuctuation (T=0) and consists of minimiz- 
ing the energy of one tetrahedron with respect to the 
spins coordinates (the eight angles defining the directions 
or the 4 spins of one tetrahedron). 

Doing so, one finds that there is a continuous degener- 
acy of states which minimize the energy of one tetrahe- 
dron. These states can be classified in two sets. The first 
one is made of the coplanar states obtained in the monte 
carlo simulations and have a continuous global degree of 
freedom which is a global rotation in the plane (rota- 
tion around z on figure ^ . The second set of lowest en- 
ergy states contains non coplanar states which can be de- 



scribed starting from a coplanar structure. Starting from 
the coplanar structure of figure |S1 one can parametrize 
the non coplanar states as follows : 

{C0S6'C0S ((^ — f ) f COS0COS ( — (y5 + f ) 

cos6'sin(v3- f) S2 = < cos6lsin(-¥3+ f) 
sin (61) [ - sin (0) 

{COS0COS (-(p - ^) ( COS 9 COS (if + ^) 

cosesm{-ip-^) S4=< cos6lsin((^+ ^) 
- sin {9) [ sin (9) 

where the spins are labeled as in figure |21 and where ip 
and 9 are not independant : 

6* = arctan(\/2siniy9) (2) 

As soon as Q is true, the corresponding state is one of 
the degenerate ground states. The state represented on 
figure ini corresponds to 9 = if = 0. In ^ 9 is restricted 
to [— f , f], however, starting from a state equivalent to 
the one of figure but where the spins are coplanar in 
the (zx) or (yz) plane, one can write down the same kind 
of parametrization of lowest energy states. Note that 
the degeneracy due to the free choice of one of the two 
angles {9 or ip in equation Q) corresponds to a global 
degree of freedom for the pyrochlore lattice, so that the 
macroscopic degeneracy of the pyrochlore antiferromag- 
net without DMI is lifted and thus a phase transition due 
to DMI is expected already on a mean field level. 

The conclusion of the mean field treatment is thus that 
the coplanar state represented on figure |S1 is a ground 
state as well as the equivalent coplanar states in the (zx) 
and (yz) planes, and these three states can be obtained 
one from another by distorting continuously the magnetic 
structure while staying at the minimum of the energy. 
However, the intermediate states are not coplanar. The 
structure (in phase space) of the lowest energy states is 
sketched on figure [7| As we shall see, the planar and non 
coplanar states are not equivalent as soon as temperature 
is not zero : they all minimize the energy but the thermal 
fluctuations will favour the planar states. 

Finally, the magnetic structure for the direct and indi- 
rect cases are very different. There is however no reason 
for them to be related on a frustrated (non-bipartite) lat- 
tice. For instance the magnetic strucutres for J > and 
J < are very different on the pyrochlore lattice, the fer- 
romagnetic system being magnetically ordered, whereas 
the antiferromagnetic has a spin-liquid ground state. 



3. Order by disorder 

This discrepancy between the two approaches (mean 
field and monte carlo) is interpreted as an entropic ef- 
fect : the mean field approach neglects the thermal fiuc- 
tuations and identifies the ground states with the minima 



FIG. 6: Ground state in the case of indirect DMI. The ground 
state for the whole pyrochlore lattice is a q = structure so 
that only one tetrahedron is represented. Similar structures 
in the zx and yz planes are degenerate. Other non-coplanar 
states have the same energy but do not participate in the low 
temperature properties (order by disorder, see text). 



of the energy. However, the fluctuations around the dif- 
ferent minima are not equivalent. Some of them are en- 
tropically favorable, and since at finite temperature the 
system minimizes its free energy, the system will prefer- 
ably fluctuate around the minima where the entropy is 
the higher, thus introducing a difference between the de- 
generate ground states at T = 0. Those states found in 
monte carlo simulations are a subset of the one found by 
energy minimization. However, this "order by disorder" 
is not responsible for the ordering of the system as it 
only selects coplanar states by removing a global degree 
of freedom (associated with equation ||2Jl) and fixes the 
relative orientation of the four Neel sublattices, but does 
not remove a macroscopic number of degrees of freedom. 

In order to be more quantitative, a coplanarity pa- 
rameter (which measures to what extent the system is 
coplanar) was defined as follows : 



FIG. 7: Schematic view of the low lying states in the presence 
of indirect DMI. (xy), (zx) and (yz) are coplanar states. En- 
ergetically, these states are equivalent to non coplanar states 
and the magnetic structure can be continuously distorted to 
explore all these energy minima, as represented by the lines 
of the picture. Starting from point A, the thin line can be 
described by varying 8 and (p according to Eq. Q (point A : 
9 ^ tp = 0; point B : 9 = (p = —n/4; point C : 6 = (p = 7r/4). 
At finite temperature, the coplanar states are selected by the 
thermal fluctuations which is represented schematically by 
broad lines. 
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FIG. 8: Evolution of the coplanarity parameter C (see equa- 
tion Q) with temperature. C = 1 corresponds to states where 
all the spins lie in one of the (xy), (zx) or (yz) planes. These 
coplanar states are selected at low temperature. 



C = 1 in any of the three coplanar states ((xy), (zx) 
or (yz)), and C = in a paramagnetic phase. Figure 
|S1 represents the evolution of the coplanarity parameter 
with temperature. One clearly sees that the coplanar 
states are selected. From the monte carlo simulations 
done so far, all the coplanar states (obtained by a rotation 
around the z axis of figure do no seem to be exactly 
equivalent, however we could not find a clear evidence of 
a (possible) further order by disorder which would select 
one particular direction in the selected planes. 



4.. Pyrochlore compounds 

The magnetic structure which appears in the pres- 
ence of indirect DMI was not observed experimentally as 
clearly as in the direct case, however, there are two com- 
pounds whose magnetic structure shares a lot of proper- 
ties with the one obtained here. 

One of these compounds is very similar to FeFa in 
which the direct DMI structure is observed. This com- 
pound is NH4Fe^^Fe^^^F6 and differs from FeFs by the 
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presence of NHj which imphes both Fe^+ and Fe^"*" mag- 
netic ions with the same concentration. The latter has 
one extra electron compared to the half filling case of 
Fe'^"'". The observed magnetic structure is similar to the 
one of figure IHl except that the 90° angle between the 
spins on the figure is ^76° in the compound [2^. Since 
the Fe^+ ion is no more isotropic, it is not surprising 
that the observed magnetic structure is a little bit dif- 
ferent from the perfect DMI structure since one expects 
other sources of anisotropics (like single ion anisotropics) 
in such a system. This however suggests that DMI may 
also be present in NH4Fe^^Fe^^^F6 but with a change of 
sign of the DMI compared to FeF^. 

The other compound whose low temperature long 
range ordered structure may share similarities with the 
one found with indirect DMI is the spinel compound 
ZnCr2 04, although the experimental structure is not 
uniquely determined |29j |. One of the candidates struc- 
tures has a planar anisotropy in the {001} planes and the 
magnetic structure on one tetrahedron has the spins anti- 
parallel two by two and otherwise perpendicular to each 
other as shown in figure |B| The experimentally observed 
wave vector is however not q = 0. The discrepancy could 
be due to the structural distortion which occurs in this 
compound and which lowers the symmetry. This is not 
considered in the present work. 



C. CU4O3 (paramelaconite) 

CU4O3 is a mineral compound in which the spin-i 
Cu^"*" ions form a tetragonally distorted pyrochlore lat- 
tice. The crystallographic structure is elongated along 
the (001) direction of the cubic cell. On the contrary of 
what is predicted for the (undistorted) pyrochlore struc- 
ture with nearest neighbour antiferromagnetic interac- 
tions, this compound orders magnetically at low tem- 
perature 24]. On-site anisotropics should be small due 
to the isotropy of the spin-i and can not give rise to a 
critical temperature of 40 K. It seemed natural to try 
and explain the ordering by introducing different cou- 
pling constants depending on the direction of the bond, 
due to the distortion, as well as small interactions be- 
tween some of the second nearest neighbour which should 
be present in this compound ^2^. However, the experi- 
mentally observed magnetic structures j24| could not be 
reproduced, not even the correct wave vector. Exper- 
imentally, the magnetic structure is not fully resolved 
and there are two possible candidates Introducing 
DMI in accordance with the crystallographic structure, 
together with the previously tried isotropic exchanges, 
we were able to reproduce one of the candidates for the 
magnetic structure. More details about this analysis were 
published elsewhere [30l |. 



D. Quantum fluctuations 

All the work presented in this paper treats the mag- 
netic moments as classical 3D vector spins of fixed length. 
This is usually only a good approximation for high spin 
quantum numbers but is believed to produce at least 
qualitatively good results for the systems under study 
here, and we now give some clues which support our ap- 
proach. First it was shown in a similar study on the 
kagome lattice that quantum effects are quite small in 
the presence of DMI interactions [sij. Second, the py- 
rochlore lattice is three dimensional and so the quantum 
effects are expected to be even smaller than in the kagome 
lattice (2D). Third, if the ordering of the paramelaconite 
(CU4O3) is indeed due to DMI as is proposed, then it is an 
experimental fact that quantum fluctuations are not very 
important in such systems since the magnetic moment of 
the S = ^ Cu^+ ion is still significant (2j|. Finally, a 
related work |l6j| considered the extreme quantum case 
of spins S = ^ (using other approximations) and both 
approaches give similar (although not exactly identical) 
results. 

The classical treatment presented here should thus be 
relevant at least qualitatively even for the extreme quan- 
tum case of spins S — However, quantum fluctuations 
were shown to renormalize strongly the critical temper- 
ature in a 2D frustrated system of spins S = ^'sS], and 
the critical temperature should also be renormalized by 
quantum fluctuations in the case studied here, although 
this effect should be weaker because the system is 3D. 



IV. CONCLUSION 

In this paper, the highly frustrated pyrochlore antifer- 
romagnct with additional DMI is studied. The relevance 
of such interactions for this system is shown and the ge- 
ometry of the D vectors which define them is deduced 
from symmetry arguments. Two kinds of DMI can exist 
in a pyrochlore lattice and the consequences of both on 
the magnetic properties are investigated. They are found 
to be in deep contrast with those of the nearest neigh- 
bour Heisenberg antiferromagnet on this lattice which 
has a spin disordered state at low temperature. Indeed, 
as soon as DMI are present, a phase transition to a mag- 
netically long range ordered state will take place at a 
temperature which is of the order of magnitude of the 
DMI. The exact magnetic structures are described and 
the relevance of this ordering scenario is discussed for 
compounds in which these structures are experimentally 
observed. 
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